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An empirical analysis of interest rates in money and capital markets is performed. 
We investigate a set of 34 different weekly interest rate time series during a time 
period of 16 years between 1982 and 1997. Our study is focused on the collective 
behavior of the stochastic fluctuations of these time-series which is investigated by 
using a clustering linkage procedure. Without any a priori assumption, we individu- 
ate a meaningful separation in 6 main clusters organized in a hierarchical structure. 
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1 Introduction 

Since long time financial data have been widely studied by economists, math- 
ematicians and, more recently, by physicists [1,2,3]. The variations of these 
financial time series can be seen as stochastic processes where a set of fi- 
nancial quantities is varying in time as consequence of underlying economic 
changes. The present availability of enormous sets of financial data allows to 
get increasingly important insights on the complex behavior of these systems 
starting from empirical studies. These investigations are leading to more and 
more accurate results on risk assessment and search for market imperfections. 
One of the important points in these analyses is to individuate similarities 
and specificities among the analyzed financial time series. This search has 
been widely exploited for stocks price changes whereas interest rates have 
been less investigated [4,5,6,7,8,9]. For several economic reasons, interest rates 
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Table 1 

Interest rates and standard deviations in the time period 1982-1997. 
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and bonds have very similar statistical behavior in time or, in other words, 
they are all highly correlated . Their multivariate dynamics has been studied 
with a correlation-based clustering procedure in a set of US treasury secu- 
rities where an underlying hierarchical structure has been detected [7]. Here 
we investigated a partially different and less homogeneous set to evaluate the 
degree of hierarchal organization among different time series in a diversified 
group of bonds. 



2 An empirical analysis on interest rates 



We investigate weekly data for 34 selected interest rate time series recorded 
in the Federal Reserve (FR) Statistical Release database [10]. In the following 
we will indicate these time series with the symbol where t is the current 

date and i is a number which labels the different time series (see Tab.l). The 
different interest rate time series analyzed are: The Federal funds rate (FED); 
State & local bonds (SLB); Commercial Paper (CP); Finance Paper placed 
directly (FP); Bankers acceptances (BA); The rate on certificates of deposit 
(CD); (Note that in these cases the numbers 1, 3 and 6 stand for maturity 
dates of 1, 3 and 6 months.); The yields on Treasury securities at 'constant 
maturity' (TC) (in particular the TC at 3 and 6 months (TC3M, TC6M) and 
1, 2, 3, 5, 7, 10, and 30 years (TC1Y-TC30Y) maturities; The Treasury bill 
rates (TBA) with maturities of 3 and 6 months (TBA3M, TBA6M, TBS3M, 
TBS6M) and 1 year (TBSIY); The Treasury long-term bond yield (TClOP); 
The Eurodollar interbank interest rates (ED) with maturity dates 1 month, 3 
months and 6 months (EDIM, ED3M, ED6M), respectively; The Corporate 
bonds Moody's seasoned rates (AAA, BAA) and The Conventional mortgages 
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Fig. 1. Left) Interest rates fi{t) as function of t for i = 1 — 34 (grey lines) and their 
average f{t) (black line). Right) Cumulative distributions for the probabilities of 
the interest rates fluctuations, (a) ^'+(x) and (b) ^-(x). (The black lines are their 
averages.) 



rates (CM). Their characteristics can be found in [10]. Unless differently stated, 
we report weekly data obtained from unweighted averages of daily data ending 
on Friday. 



3 Fluctuations 



The interest rate time series, fi{t) v.s. t are shown in Fig.l, where their av- 
erage f{t) = J2i fiif)l'i^ is also shown. It is evident from Fig.l that all these 
data follow very similar trends in time and they lay in a narrow band around 
f{t). The interest rate fluctuations are analyzed by studying the changes in 
their values from one week to the following week: A/j(t) = fi{t + At) — /j(t) 
where At = 1 week. The quantities Afi{t) show stochastic fluctuations around 
the zero with similar behaviors for all the interest rates. These fluctuations 
are analyzed in the 'tails' region by computing the cumulative distributions 
\l/±(±a;) (x > 0), a quantity which tells us the probability to find a weekly 
change which is larger than x (+), or smaller than —x {—). It is defined as: 
*+(x) = 1 - S-ooP{C)d^ and \l/_(x) = Sl^p{C)d^ with p(^) being the proba- 
bility density distributions of Afi(t) (see Fig.l (Right)). These distributions 
are highly leptokurtic and are characterized by non-Gaussian profiles. The 
standard deviation of Afi{t) is defined as: = \J rp^^rp^ Y^LT^i^fii'^) ~ {^f)Y 
where Ti and T2 delimit the range of t, and (A/) is the average over time 
of Afiit) (which tends to zero for T2 — Ti ^ 00). We compute the standard 
deviations of A/j for each interest rate series for the whole period 1982-1997 
(see Tab.l) and for each year (see Fig.2 (Left)). We can observe an overall 
decreasing trend of cxj in the time period 1982-1997 with similar fiuctuations 
for all the interest rates series except for FED. 
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Fig. 2. Left) Standard deviations of /S.fi{t) for all the rates analyzed (i = 1..34) as 
function of the years between 1982-1997. Right) Hierarchical tree obtained from the 
correlation coefficients of the 34 interest rates fluctuations time series A/j(t) in the 
time period 1982-1997. (On the x-axis are reported the i values and on the y-axis 
the ultra- metric distances.) 

4 Cluster analysis and discussion 



To understand the geometrical and topological structure of the correlation co- 
efficients, we use the metric distance dij between the series A/j and Afj which 



is defined in [11] and used for financial time series in [12]: dij = ^2{1 — Cij) 
with Cij the correlations among the i,j interest rates weekly changes: 

_ (A/.A/,) - (A/.) (A/,) 

(TiCTj 



where the symbol (...) denotes a time average performed over the investigated 
time period. The correlation coefficients are computed between all the pairs 
of indices labeling our interest series. Therefore we have a 34 x 34 symmetric 
matrix with = 1 on the diagonal. By definition, Qj- is equal to zero if the 
interest rates series i and j are totally uncorrected, whereas = ±1 in the 
case of perfect correlation/anti-correlation. Therefore c/jj can vary between 
to 2. We determine an ultra-metric distance di^j which satisfies the first two 
properties of the metric distance and replaces the triangular inequality with 
the stronger condition: d^j < max[(ij^fc, rf^ j], called 'ultra-metric inequality'. 
Once the metric distance dij is used, one can introduce several ultra-metric 
distances. Mantegna et. al have used the 'subdominant ultra-metric', obtained 
by calculating the minimum spanning tree connecting several financial time 
series [13,14,15,16]. Here we consider a different ultra-metric space that em- 
phasizes the cluster-structure of the data. In our case, a "cluster" is a set of 
elements with relative distances d^ j which are smaller than a given thresh- 
old distance 5, whereas disjoined clusters have some elements which are at 
distances larger than 6. We define the ultra-metric distance dij between two 
distinct elements i, j belonging to two different clusters as the maximum metric 
distance between all the couples of elements in the two clusters [8]. The link- 
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Fig. 3. Left) Ultra-metric distances 5n and 6p at which the clustering process begins 
and ends, as function of the years between 1982-1997. Right) Cluster-structure per- 
sistence in the period 1982-1997. In the first column the interest rates are indicated. 
In the second column, the grey tones distinguish the different clusters as resulting 
from the analysis over the whole time period. The other columns refer to the cluster 
analysis over the four selected time periods, namely 82-85, 86-89, 90-93, 94-97. The 
numbers inside each cluster refer to the ultra-metric distances at which each cluster 
starts and ends its clusterization. 



age procedure yields to a hierarchical graph as shown in Fig. 2 (Right), which 
refers to the 34 time series in the whole period 1982-1997. The clustering pro- 
cess starts with a nucleation between TCIY and TBSIY at the ultra-metric 
distance 6n = 0.087. The clustering ends when all the series merge in a unique 
large cluster at the ultra-metric distance Sp = 1.3083. At this distance all the 
interest rates with maturity dates smaller or equal than 6 months (already 
merged with the FED at d = 1.18) make a single cluster with another cluster 
composed of interest rates with maturity dates larger or equal than 1 year. 
The same analysis performed on each year, gives comparable 5.„ and 6p values 
which are plotted in Fig. 3 (Left). Let now consider the intermediate region by 
analyzing the cluster evolution at the threshold distance 5 = l/\/2 = 0.707..., 
which is half the way between completely uncorrelated series (cjj = and 
dij = v^) and completely correlated ones (qj = 1 and dij = 0). At this 
threshold distance, the cluster analysis on the whole data set (1982-1997) 
leads to 6 clusters and 3 isolated elements, as one can see from Fig. 2 (Right). 

The corresponding interest rates associated with these clusters are summarized 
in first column (1982-1997) of Fig.3 (Right) where the ultra-metric distances at 
which the nucleation process starts and ends for each cluster are also indicated. 
As one can see, the empirical analysis allow us to distinguish several different 
clusters that gather together meaningful quantities: 

- all the interest rates with maturities equal to 1 months; 
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- all the interest rates with maturities 3 and 6 months with distinctions for the 
Treasury securities at 'constant maturity' (TC), Treasury bill rates (TBA) 
and Treasury bill secondary market rates (TBS); 

- all the interest rates with maturities between 1 and 3 years; 

- all the interest rates with maturities larger than 3 years. 

Fig. 4 reports the plot for the interest rates time series fi{t) grouped into the 
different sets retrieved from the cluster analysis described above. For some 
of the series the data-collapse is impressive, indicating that the correlations 
inside the clusters are strong in any part of the analyzed period. It is therefore 
interesting to investigate whether a cluster structure, similar to the one ob- 
tained for the period 1982-1997, could be retrieved from an analysis on shorter 
periods. This is of course a delicate point since the fragmentation of the data 
sets will increase the fluctuations due to the noise. We choose to divide the 
whole period 1982-1997 in four smaller periods of 4 years. The results are re- 
ported in Fig. 3 (Right) where it is evident how the cluster structure is mostly 
conserved (and partially modified) in this 4-years period analysis. In conclu- 
sion, from the analysis of different kinds of interest rates in money and capital 
markets, referring to government, private, industries securities and commit- 
ments, we have shown how the used clustering linkage procedure is useful to 
detect differences and analogies among these tangled correlated data. 
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